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1 INTRODUCTION 


ABSTRACT 

A recent survey of the inner 0.35 x 0.35 pc of the NGC 2024 star forming region 
revealed two distinct millimetre continuum disc populations that appear to be spatially 
segregated by the boundary of a dense cloud. The eastern (and more embedded) 
population is ~ 0.2 — 0.5 Myr old, with an ALMA mm continuum disc detection rate 
of about 45 per cent. However this drops to only ~ 15 per cent in the 1 Myr western 
population. When presenting this result, van Terwisga et al. (2020) suggested that the 
two main UV sources, IRS 1 (a B0.5V star in the western region) and IRS 2b (an O8V 
star in the eastern region, but embedded) have both been evaporating the discs in the 
depleted western population. 

In this paper we report the firm discovery in archival HST data of 4 proplyds and 4 
further candidate proplyds in NGC 2024, confirming that external photoevaporation 
of discs is occurring. However, the locations of these proplyds changes the picture. 
Only three of them are in the depleted western population and their evaporation is 
dominated by IRS 1, with no obvious impact from IRS 2b. The other 5 proplyds 
are in the younger eastern region and being evaporated by IRS 2b. We propose that 
both populations are subject to significant external photoevaporation, which happens 
throughout the region wherever discs are not sufficiently shielded by the interstellar 
medium. The external photoevaporation and severe depletion of mm grains in the 
0.2-0.5 Myr eastern part of NGC 2024 would be in competition even with very early 
planet formation. 
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2015; Rodriguez et al. 2018; Winter et al. 2018b; Concha- 
Ramírez et al. 2019a; Moore et al. 2020), chemical enrich- 


Planet formation is increasingly being considered in the 
wider context of star forming regions (for a recent review 
see Parker 2020). The environmental impact on discs and 
planets from the surroundings may take place through ex- 
ternal photoevaporation (e.g. Stórzer & Hollenbach 1999; 
Adams et al. 2004; Facchini et al. 2016; Guarcello et al. 2016; 
Kim et al. 2016; Ansdell et al. 2017; Eisner et al. 2018; Ha- 
worth & Clarke 2019; Nicholson et al. 2019; Reiter & Parker 
2019; Concha-Ramírez et al. 2019b; Winter et al. 2020a), 
gravitational encounters (e.g. Cabrit et al. 2006; Dai et al. 
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ment, such as through short lived radionuclides from super- 
novae (Kuffmeier et al. 2016; Close & Pittard 2017; Lichten- 
berg et al. 2019) and cosmic ray ionisation (Kuffmeier et al. 
2020). In terms of direct impact on the physical character- 
istics (mass/radius) of an established disc, external photo- 
evaporation is generally thought to dominate over encoun- 
ters (Scally & Clarke 2001; Winter et al. 2018a). Cosmic 
rays may affect the initial conditions of discs through the 
strength of magnetic braking (Kuffmeier et al. 2020). 


We have growing statistical evidence for the impact of 
radiation environment on discs (e.g. Fang et al. 2012; Mann 
et al. 2014; Maucé et al. 2016; Guarcello et al. 2016; Ans- 
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dell et al. 2017; Eisner et al. 2018). However an ongoing 
issue is that we have very limited sample of individual discs 
with observed photoevaporative outflows. T'he most intense 
study has been of photoevaporating individual discs in the 
Orion Nebular Cluster (ONC) - the famous “proplyds” with 
cometary morphologies consisting of a cusp and tail, as il- 
lustrated schematically in Figure 1 (McCaughrean & O'dell 
1996; O’dell 1998; Henney & O’Dell 1999; Bally et al. 2000; 
Henney et al. 2002). Other photoevaporating discs are in 
similarly strong, or stronger, UV environments including in 
Pismis 24 (5 candidate proplyds in a UV environment even 
stronger than the ONC - Fang et al. 2012) and Carina (e.g. 
Smith et al. 2010) 


Recently Kim et al. (2016) discovered 7 proplyds in 
NGC 1977, a UV environment around 2 orders of magni- 
tude weaker than the ONC. This provided important first 
direct evidence for external photoevaporation operating in 
weaker UV environments. However, more examples are re- 
quired across a larger range of UV environments to truly 
understand the photoevaporation process in action. To this 
end we require both more proplyd detections and to under- 
stand evaporation signatures so that we can detect a wind 
when it is more subtle (e.g. in CI and CO lines - Haworth 
& Owen 2020). 


'The stellar cluster NGC 2024 was recently observed 
with ALMA by van Terwisga et al. (2020), who found that 
there were two disc populations that are similar in age, but 
separated both spatially and in terms of disc fraction. The 
slightly younger (~ 0.2— 0.5 Myr) population, within and on 
the eastern side of a dense molecular cloud, has an ALMA 
1.3mm disc fraction of 45 + 7 per cent (the rms noise of 
observations was 0.058 mJy beam ! and lowest “detection” 
dust mass was 2.5 + 0.70 Ma van Terwisga et al. 2020) and 
masses similar to 1-3 Myr old star forming regions such as 
Lupus. The slightly older (~ 1 Myr) population lies on the 
western side of the cloud and has an ALMA mm disc fraction 
of only 15+ 4 per cent. Being less shielded and containing 
the BO star IRS 1 it was speculated that external photoevap- 
oration may have been responsible for depleting the western 
disc population. However it is also possible that discs in the 
western region simply began with lower masses. Note that 
the near infrared disc fraction over a much larger extent of 
the region (6.7 x 6.7’ vs 2.9 x 2.9") is higher, at around 85 per 
cent (Haisch et al. 2001a,b) 


Kuffmeier et al. (2020) recently suggested that higher 
cosmic ray ionisation rates can enhance the impact of 
magnetic braking and lead to initially smaller (and hence 
possibly shorter lived) discs. They suggested that this 
might provide an alternative explanation for the different 
sub-mm detection fractions in NGC 2024. It is important 
to understand what has had such a large impact on these 
two disc populations, both to understand this region and 
the planet formation process more generally. Determining 
the impact of external photoevaporation in this type of 
environment is therefore key. 


Motivated by van Terwisga et al. (2020), here we search 
archival Hubble Space Telescope (HST) data for photoevap- 
orating discs in NGC 2024. We discover 4 proplyds and 4 
further candidates which demonstrate that external photo- 
evaporation is indeed operating throughout the region. 


Anatomy of a proplyd 
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Figure 1. A schematic of a cometary proplyd. The disc is being 
irradiated by a UV field represented by the parallel arrows. This 
drives material from the disc into a bright cusp on the side facing 
the UV source and a tail on the side away from the source. When 
measuring the size of the proplyd (the ionisation front radius in 
our case) we use the radius of a circle that traces the cusp. 


2 OVERVIEW OF UV SOURCES AND STRUCTURE 
OF NGC 2024 


A mosaic overview of NGC 2024 in the HST F160W filter 
(~ 1.6um) constructed with data from the MAST" science 
archive (PI: Arce) is given in Figure 2. Overlaid are Her- 
schel PACS blue 70um contours from 1-6 Jy pixel ! using 
data from the ESA Herschel Science Archive? (PI: Megeath). 
'The diagonal line at the western edge of the cloud separates 
the younger (~ 0.2 — 0.5 Myr) population from the 1 Myr 
old population that exhibit relatively depleted mm-fluxes 
(Haisch et al. 2000; Getman et al. 2014; van Terwisga et al. 
2020). The star symbols denote the positions of the main 
UV sources: IRS1 (B0.5V star), IRS2b (08V star) and IRS 
2 (B star), which are summarised in Table 1. The circles 
denote the positions of proplyds introduced in this paper 
(discussion of these is left for sections 5 and 6. 

NGC 2024 is the youngest cluster in Orion and rich- 
est in star formation in Orion B (Meyer 1996; Levine et al. 
2006). There are two main UV sources. IRS 1 (a B0.5V star 
— Garrison 1968) is not embedded and resides in the ~ 1 Myr 
region which van Terwisga et al. (2020) conclude has a de- 
pleted mm disc population. Given our interest in proplyds it 
is worth noting that IRS 1 is similar to 42 Ori in NGC 1977, 
where 8 proplyds have now been discovered (Bally et al. 
2012; Kim et al. 2016). 

IRS 2b (an O8V star Bik et al. 2003) is in the eastern 
(~ 0.2— 0.5 Myr) population (van Terwisga et al. 2020). IRS 
2 is an embedded B-type star in close projected proximity 
to IRS 2b and is expected to be of secondary importance in 
terms of its contribution to the UV flux in the region. 

For some time only IRS 1 and 2 were known in NGC 
2024, giving rise to the problem that neither have hard 
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Figure 2. An overview of NGC 2024. The main image is a 1.6 4m F160W HST mosaic taken from MAST (Arce et al. 2016). Overlaid are 
Herschel PACS blue 70um contours, from 1-6 Jy/pixel (Stutz et al. 2013). The 60" grid is 0.12 pc at a distance of 414 pc (Menten et al. 
2007; Bailer-Jones et al. 2018). Star symbols denote the main UV sources (IRS1 and IRS2). Circles denote the three proplyds introduced 
in this paper. The diagonal dashed line separates the two populations of marginally different age, but substantially different disc fraction 


(van Terwisga et al. 2020). 


enough a spectrum to explain the degree of hydrogen and he- 
lium ionisation (Bik et al. 2003). The discovery of the more 
massive IRS 2b appears to resolve this (e.g. Thompson et al. 
1981; Kruegel et al. 1982; Barnes et al. 1989; Bik et al. 2003; 
Kandori et al. 2007; Burgh et al. 2012), though it is still not 
ruled out that there are perhaps other UV sources on the far 
side of the dense cloud. Kruegel et al. (1982) and Bik et al. 
(2003) respectively suggest ~ 104% and 7.31.2 x 10^ ionis- 
ing photons per second are being absorbed by the gas to give 
the observed radio continuum signal. Note that these values 
are a lower limit since dust absorption of ionising photons, 
or ionising photons escaping from the region, will not ionise 
the gas and contribute to the radio continuum. IRS 1 can- 
not explain all of the ionising flux expected from the radio 
continuum, but the exact value of the relative contributions 
from IRS 1/ IRS 2b is unknown. 


MNRAS 000, 000-000 (0000) 


In many star forming regions a single ionising source is 
dominant. However, NGC 2024 has a complicated structure 
with a lot of gas and dust in the region, so it isn't obvi- 
ous that IRS 2b necessarily dominates the UV of the entire 
region. 


'The cusp of proplyds point towards the UV source re- 
sponsible for evaporating them (see Figure 1), so the de- 
tection of proplyds in NGC 2024 would also shed light on 
where the different UV sources dominate over one another. 
Furthermore, if proplyds were discovered pointing towards 
none of the known UV sources, it would give clues to the 
location of additional massive star(s). 
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Table 1. The main known UV sources in NGC 2024. 


UV RA DEC Spectral 
Source class 


IRS 1 05 41 37.8540 
IRS 2 05 41 45.81 
IRS 2b 05 41 45.50 


-01 54 36.6256 B0.5V 
-01 54 29.8 B 
-01 54 28.7 O8V 


3 THEORETICAL POTENTIAL FOR PROPLYDS 


For proplyds to be visible and resolvable requires a com- 
bination of sufficient EUV irradiation and a high-enough 
mass-loss rate. The radius Rire of the ionisation front can be 
approximated by assuming a balance between EUV flux and 
recombinations in the flow (see Johnstone et al. 1998): 


"m 2/3 
10-8 Ms —) 


N —1/3 d 2/3 
x (usi 2) (5) au (1) 
where Nr, is the EUV photon count from the (external) 
ionising source, d is the distance to that source, and Mind 
is the wind driven mass-loss rate from the disc. 

We estimate the radiation field at any given point in 
NGC 2024 as a function of the distance (in practice, the 
projected separation) from the two known main UV sources: 
IRS 1 and IRS 2b (see section 2/Table 1). 

The intrinsic hydrogen ionising photon count from 
these sources, following Sternberg et al. (2003), is 
logio(Nry msi) = 47.71 and logo (Nzy,ms2b) = 48.75 for 
IRS 1 and IRS 2b respectively. This does give a higher in- 
put ionising flux than inferred from the radio continuum 
(e.g. Bik et al. 2003), but as discussed in section 2 the radio 
continuum does only provide a lower limit. In the discussion 
section 6.1 we will also consider values a factor 5 lower for 
the ionising photon emission rate from IRS1/IRS2b. 

As a function of distance d (projected separation) we 
hence have an ionising photon flux of 


Rig Evi 1200 ( 
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for IRS1 and IRS2b respectively. 

The ionising flux is important for estimating mass loss 
rates as per equation 1. It is also useful to estimate the 
magnitude of the FUV radiation, which we do to first or- 
der by integrating over a blackbody spectrum, assuming 
stellar temperature and luminosity from Sternberg et al. 
(2003). This gives Lruv rsı ^ 1.4 x 10?? erg s^! and 
Lruv,ns2b ~ 3.4 X 1038 erg s^! for IRS1 and IRS 2b re- 
spectively. IRS1 is in a relatively open part of the region, so 
there is comparatively little extinction to the stellar neigh- 
bours in the cluster. Extinction may play a large role in the 
vicinity of IRS2b though, with many of the nearby stars be- 
ing embedded. Assuming the no-extinction case provides an 
upper limit, the FUV radiation field strength as a function 


of distance from the UV sources is 


d -2 
F < 7.34 x 10* | ——— 4 
UVirs1 SS 7.34 x 10 (st) Go ( ) 
and 

d —2 
F 1.8 x 10 | — 
UVirsap < 1.8 x 10 (st) Go (5) 


for IRS1 and IRS2b respectively. where Go is the Habing 
unit of the radiation field’. 

Given that the region studied by van Terwisga et al. 
(2020) and illustrated in Figure 2 is around half a parsec in 
size, these ionising and FUV flux estimates are high enough 
that we theoretically expected them to drive significant pho- 
toevaproative winds from discs in the region, as long as they 
are not deeply embedded (e.g. Johnstone et al. 1998; Ha- 
worth et al. 2018b). 


4 ARCHIVAL HST OBSERVATIONS 


We searched for historical Hubble Space Telescope (HST) 
observations towards NGC 2024 in the MAST* science 
archive and Hubble Legacy Archive? (HLA). 


41 WFC2 Ha 


We studied Wide Field Camera 2 Ha images from the HLA. 
For these all we required was the morphological identifica- 
tion of proplyds, their coordinates and measurements of the 
ionisation front size. We therefore used the data directly as it 
appears in the science archive, with no further manipulation. 
We use a single field that is image 1 from Proposal ID 5983, 
F656N, PI: Stapelfeldt 1995. The pixel size is 0.1/’which cor- 
responds to 41 au at a distance of 414 pc (Menten et al. 2007; 
Bailer-Jones et al. 2018). 


4.2 NICMOS Paschen o 


The main dataset that we use was originally obtained in 2003 
by Stapelfeldt & Hines (Proposal ID: 9424, PI Stapelfeldt). 
These used camera 2 of NICMOS in filters centered on 
Paschen a (F187N) and the nearby (in wavelength) contin- 
uum (F190N) to image fields centered on 16 compact VLA 
sources in the region. Subtracting the F190N image from 
the F187N therefore leaves continuum subtracted Paschen a 
emission. The pixel size is 0.075" which corresponds to 31 au 
at a distance of 414 pc. 

A 3-position spiral dither pattern was executed for each 
object, with 5” offsets to mitigate against bad pixels, resid- 
ual cosmic rays and other imaging artifacts within 10” of the 
centre of the field of view. Three positions (the usual rec- 
ommended pattern uses four positions) were used to maxi- 
mize the serendipitous field of view and yet still fit within 
the allocated observing time. After processing through the 
STScI pipeline, the three dithered images were shifted into 


3 Note that an alternative measure of the UV field is the Draine 
X, which is 1.71 times larger than the Habing unit. Hence 1 Go is 
~ 0.585 x 
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common world coordinates and median combined. The re- 
sulting mosaic images were used to construct Paschen a line 
images by subtracting the F190N image scaled by the ap- 
propriate bandpass (the stellar photosphere is assumed to 
be flat across the F187N and F190N bandpasses). 


5 PROPLYDS IN NGC 2024 
5.1 Overview 


We have firmly identified 4 proplyds and have 4 other can- 
didate proplyds in NGC 2024. Two of these were found in 
the F656N filter (Ha) and the other 6 in F187N-F190N 
(Paschen o). Three proplyds are in the western region (the 
^ 1Myr population with only ~ 15 per cent of discs de- 
tected in the mm-continuum by van Terwisga et al. 2020) 
and pointing towards the B0.5V star IRS 1. The other 5 are 
in the ~ 0.2— 0.5 Myr higher mm disc fraction region, point- 
ing towards the O8 star IRS2b. These will all be discussed in 
detail below, but for quick reference, image galleries of the 
proplyds are given in Figures 3-6 and a summary of their 
main parameters is given in Table 2. The locations of the 8 
propylds and the main UV sources are also provided in the 
overview of the region in Figure 2. 

We next introduce the basic properties (morphology, 
location, environment) of each proplyd in turn and in the 
discussion turn our attention to the region as a whole. 

The ionisation front radii were estimated using the ra- 
dius of a circle manually drawn over the cusp of the proplyd 
(see Figure 1). Uncertainties are based on circles that trace 
the inner and outer edge of the cusp. 


5.2 Proplyd 1 


Proplyd 1, discovered in F656N H a (and preliminarily in- 
troduced in the proceedings of Stapelfeldt et al. 1997) is 
shown in Figure 3. It is at a projected separation (assuming 
a distance of 414 pc, Menten et al. 2007; Bailer-Jones et al. 
2018) of 0.055 pc from the B0.5V star IRS 1, towards which 
its bright cusp is directed. The morphology of Proplyd 1 is 
not simply cometary, with additional Ha emission extend- 
ing to the south west. We interpret this as being from a 
jet. 

From equations 2 and 4 the UV and ionising flux is 
2.4 x 109 Go and 1.4 x 10!? s^! respectively. From an ap- 
proximate manual measurement of the radial extent of the 
cusp we assume an ionisation front radius of 144 + 62au. 
'The expected mass loss rate from equation 1 for this pro- 
plyd is hence Lye x 107" Mo yr" !. Proplyd 1 also has a 
counterpart in van Terwisga et al. (2020), where a dust mass 
of 67.2 + 0.77 Ma is inferred. 

There is also archival F187N-F190N (Paschen a) data 
for proplyd 1, which is shown in Figure 4 with Ha con- 
tours overlaid. Since the point spread function of the F187N 
and F190N images are slightly different sizes, unresolved 
Paschen a emission associated with the star has led to arti- 
facts which makes the cusp and jet difficult to identify. How- 
ever, larger scale emission from the jet and possibly also a 
more extended component of the tail are discernible. 
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5.3 Proplyd 2 


Proplyd 2 was discovered in F187N-F190N (Paschen o) emis- 
sion. The F187N, F190N (continuum) and F187N-F190N im- 
ages for Proplyd 2 are shown in the upper row of panels in 
Figure 5. Proplyd 2 is 0.037 pc from IRS 1 and an arrow is 
included on each panel pointing towards it. Like Proplyd 1, 
there are bipolar streaks of emission in addition to the typi- 
cal proplyd cusp which as in the case of proplyd 1 we again 
interpret as being due to a jet. 

From equations 2 and 4 the UV and ionising flux is 
5.4 x 109 Go and 3.1 x 10!? s^! respectively. From an ap- 
proximate manual measurement of the radial extent of the 
cusp we assume an ionisation front radius of 2063-31 au. The 
expected mass loss rate from equation 1 for this proplyd is 
hence 4.3 + 1.0 x 1077 Mo yr !. 

Proplyd 2 also has a counterpart in van Terwisga et al. 
(2020), where a dust mass of 14.6 + 0.72 Mẹ is inferred. 


5.4 Proplyd 3 


Proplyd 3 was discovered in F187N-F190N (Paschen o) emis- 
sion. The F187N, F190N (continuum) and F187N-F190N im- 
ages for Proplyd 3 are shown in the middle row of panels in 
Figure 5. Proplyd 3 is 0.057 pc from the 08V star IRS 2b, 
towards which it points. An arrow is included on each panel 
in the direction of IRS 2b. There are no additional “jet” 
features for Proplyd 3. It beautifully exhibits the classical 
cometary morphology of many ONC proplyds (e.g. O'dell 
1998). 

From equations 2 and 4 the UV and ionising flux is 5.5 x 
10° Go and 1.4 x 101? s^! respectively. From an approximate 
manual measurement of the radial extent of the cusp we 
assume an ionisation front radius of 931-31 au. The expected 
mass loss rate from equation 1 for this proplyd is hence 
2.815 x 1077 Mo yr !. 

Proplyd 3 also has a counterpart in van Terwisga et al. 
(2020), where a dust mass of 13.8 + 2.56 Mg is inferred. 


5.5  Proplyd 4 (and possibly 4b) 


Proplyd 4 was discovered in F187N-F190N (Paschen o) emis- 
sion. The F187N, F190N (continuum) and F187N-F190N im- 
ages for Proplyd 4 are shown in the bottom row of panels in 
Figure 5. Proplyd 4 points towards IRS 2b, which lies at a 
projected distance of 0.054 pc. An arrow is included on each 
panel in the direction of IRS 2b. 

Proplyd 4 exhibits a clear bright cusp with only a 
faint cometary tail. There is also an intriguing secondary 
Paschen a source which may well be an evaporating binary 
companion. The projected separation between these two ob- 
jects is only about an arcsecond, corresponding to 413 au at 
a distance of 414 pc. However given that the orbital configu- 
ration is unknown 413 au represents a lower limit on the true 
binary separation. Other binary proplyds have been discov- 
ered and, depending on the binary separation there may or 
may not be an interproplyd shell formed from the collision 
of the individual proplyd outflows (e.g. Graham et al. 2002; 
Henney 2002; Hodapp et al. 2009; Vasconcelos et al. 2010; 
Wu et al. 2013, 2018). In numerical simulations of evaporat- 
ing binary propylds, Vasconcelos et al. (2010) found that for 
separations ~ 2000 au an interproplyd shell does develop, 
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Figure 3. Proplyd 1, discovered in an archival HST WFC2 F656N (Ha) image. The star responsible for evaporating it is IRS 1, a B0.5V 
star situated only 0.055 pc away. Included in the inset box is a zoom in on the proplyd itself. The material almost perpendicular to the 


line connecting IRS1 and the bright proplyd cusp is likely a stellar jet. 
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Figure 4. The Paschen o data for proplyd 1 with the Ha con- 
tours overlaid (see Figure 3). In this case the Paschen o emission 
from the star coupled with a slightly different PSF for F187N and 
F190N leads to artifacts, but extended emission from the jet and 
possibly a longer component of the tail are still visible. 


but this is not the case for binaries of separation less than 
around 200 au. With a separation of at least 400au in the 
case of Proplyd 4, it is conceivable that there might be an 
interproplyd shell, however it is not easily discerned in the 
archival Paschen a data. 

From equations 2 and 4 the UV and ionising flux is 6.1 x 
10° Go and 1.6 x 101? s^! respectively. From an approximate 
manual measurement of the radial extent of the cusp we 
assume an ionisation front radius of 87 +37 au. The expected 
mass loss rate from equation 1 for this proplyd is hence 
2,7112 x 10 Mo yr. 

Proplyd 4 also has a counterpart in van Terwisga et al. 
(2020), where a dust mass of 26.8 + 0.87 Ma is inferred. 


5.6 Candidate Proplyd 5 


We now turn our attention to the more subtle, candidate, 
proplyds discovered in NGC2024. Although less clear, these 
all exhibit a proplyd morphology in the direction of either 
IRS 1 or IRS 2b. The candidate proplyds are also at larger 
distances, which probably explains why the proplyd mor- 
phology is less obvious. 

Candidate Proplyd 5 was discovered in F187N-F190N 
(Paschen a) emission. A F187N-F190N image of Proplyd 5 is 
shown in the upper left panel of Figure 6. Candidate Proplyd 
5 points towards IRS 2b, which lies at a projected distance 
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Table 2. Properties of proplyds in NGC 2024. Columns are the proplyd ID, coordinates, the UV source responsible for the photoevap- 
oration, projected separation from that source (assuming a distance of 414 pc), ionisation front radius, estimated mass loss rate, any 
counterpart and the corresponding dust mass from the van Terwisga et al. (2020) dataset. Note that the Paschen a pixel size is 31au at 


414 pc. 
Object RA DEC UV Projected Rit M Van Terwisga et al. Ma 
source Separation (pc) (au) (1077 Mo yr- 1) 2020 counterpart (Ma) 
Proplyd 1 05 4136.28 -01 54 24.26 IRS1 0.055 144 + 62 rey IRC215 67.2 + 0.77 
(VLA 1) (0.28 from IRS 2b) 
Proplyd 2 05 41 38.08 -01 5418.41  IRSI 0.037 206 + 31 4.31.0 IRC116 14.6 + 0.72 
(VLA 4) (0.22 from IRS2b) 
Proplyd 3 05 4145.62 -01 54 57.26  IRS2b 0.057 93 + 31 2.8*13 IRC057 13.8 + 2.56 
(VLA 20b) (0.23 from IRS1) 
Proplyd 4 05 4144.93 -01 54 54.06 IRS2b 0.054 87 + 37 2 TUE IRC59 26.8 + 0.87 
(VLA 20a) (0.22 from IRS1) 
Candidate 05 4141.53 -01 5411.5  IRS2b 0.124 Jos gate IRC 123 40.4 + 0.71 
Proplyd 5 (0.121 from IRS1) 
Candidate 05 41 42.12 -01 5405.58 IRS2b 0.11 7 +31 LISI - - 
Proplyd 6 (0.14 from IRS1) 
Candidate 05 41 42.43 -01 54 09.72 IRS2b 0.10 109 +31 PA IRC 124 42.5 + 1.79 
Proplyd 7 (0.15 from IRS1) 
Candidate 05 41 34.42 -01 53 19.49 IRS1 0.19 < 60 < 0.1 7 - 
Proplyd 8 (0.36 from IRS 2b) 


of 0.124 pc. An arrow is included in the direction of IRS 
2b. From equations 2 and 4 the UV and ionising flux is 
1.2 x 109 Go and 3.1 x 101? s^! respectively. 

Candidate Proplyd 5 is very large radius cusp of 
171438 au, which corresponds to a mass loss rate of 3.25 1 x 
107" Mo yr^! from equation 1. 

Candidate Proplyd 5 also has a counterpart in van Ter- 
wisga et al. (2020), where a dust mass of 40.4 + 0.71 Mẹ is 
inferred. 


5.7 Candidate Proplyd 6 


Candidate Proplyd 6 was discovered in F187N-F190N 
(Paschen o) emission. A F187N-F190N image of Proplyd 6 
is shown in the upper right panel of Figure 6. Candidate 
Proplyd 6 points towards IRS 2b, which lies at a projected 
distance of 0.11 pc. An arrow is included in the direction of 
IRS 2b. From equations 2 and 4 the UV and ionising flux is 
1.5 x 10° Go and 3.9 x 10? s7} respectively. 

Candidate Proplyd 6 has only a faint cusp and was pri- 
marily identified through its cometary tail. 

From an approximate manual measurement of the radial 
extent of the cusp we assume an ionisation front radius of 
77 X: 31au. This corresponds to a mass loss rate of 1.1*0 4 x 
107" Mo yr^! from equation 1. 

We could find no obvious counterpart for candidate pro- 
plyd 6 van Terwisga et al. (2020) 


5.8 Candidate Proplyd 7 


Candidate Proplyd 7 was discovered in F187N-F190N 
(Paschen o) emission. A F187N-F190N image of Proplyd 7 
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is shown in the lower left panel of Figure 6. Candidate Pro- 
plyd 7 points towards IRS 2b, which lies at a projected dis- 
tance of 0.10 pc. An arrow is included in the direction of 
IRS 2b. From equations 2 and 4 the UV and ionising flux is 
1.8 x 10° Go and 4.7 x 10? s^! respectively. 

Candidate Proplyd 7 is close to candidates 5 and 6 and 
exhibits a cusp, but no obvious cometary tail. There is a 
bright Paschen o source slightly further than candidate 7 
relative to IRS 2b. This is not a result of bad subtraction, but 
rather that this is a very bright Paschen o emitter, probably 
due to strong accretion. 

From an approximate manual measurement of the radial 
extent of the cusp we assume an ionisation front radius of 
109-31 au. This corresponds to a mass loss rate of 2.0*0$ x 
107" Mo yr^! from equation 1. 

Candidate Proplyd 7 also has a counterpart in van Ter- 
wisga et al. (2020), where a dust mass of 42.5 +1.79 M 
inferred. 


Was 


5.9 Candidate Proplyd 8 


Candidate Proplyd 8 was discovered in F656N (Ha) emis- 
sion. A Ha image of Proplyd 8 is shown in the lower right 
panel of Figure 6. Candidate Proplyd 8 points towards IRS 
1, which lies at a projected distance of 0.19 pc. An arrow is 
included in the direction of IRS 1. It is also possible that the 
tail of candidate proplyd 8 is a (very asymmetric) jet, though 
the alignment of the tail with the vector from IRS 1 is sup- 
portive of it being external photoevaporation. However, it is 
interesting to note that at least 2/3 of the proplyds in the 
older western region near IRS 1 have possible jets whereas 
none of the proplyds towards IRS 2b appear to. 
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Figure 5. A summary of the main proplyds introduced in this paper that were detected through Paschen a emission. The left hand panels 
are the Hubble F187N filter (centred on the Paschen a line) central panels are the F190N filter (the continuum) and right is F187N—F190 
which is a proxy for the continuum subtracted Paschen a line emission. In each case the arrow points towards what is believed to be the 
UV source responsible for the photoevaporation. Note that Proplyd 2 may have a bipolar outflow, which is responsible for the elongated 


emission on either side. Note that the epoch is J2000. 


From equations 2 and 4 the UV and ionising flux is 2.0x 
10? Go and 1.2 x 10!! s^! respectively. Candidate Proplyd 8 
is in the weakest UV environment of the proplyds introduced 
in this paper. 


Candidate Proplyd 8 has a very small unresolved cusp 
(it spans a maximum of about 3 pixels in diameter) and 
a bright cometary tail. We assume an upper limit on the 
ionisation front radius of 60 au. This corresponds to a mass 
loss rate of < 1 x 1078 Mo yr^! from equation 1. 


Candidate Proplyd 8 has no obvious counterpart in van 
Terwisga et al. (2020) and may be beyond the field they 
studied (it is beyond the limits of their Figure 1). 


6 DISCUSSION 
6.1 Mass loss rates and proplyd lifetimes 


Figure 7 shows the equation 1 predicted ionisation-front ra- 
dius of a proplyd as a function of distance from IRS 1 (upper 
panel) and IRS 2b (lower panel). Different lines correspond 
to mass loss rates from 10? — 10 5 Mo yr ^! using the UV 
fluxes from Sternberg et al. (2003). The shaded areas span 
a UV flux down to a factor 5 lower. Overlaid are the ioni- 
sation front radii and projected separations of our observed 
proplyds on the panel of the suspected UV source responsi- 
ble for photoevaporation (note that the estimated mass loss 
rates are all stated in section 5 / Table 2). 

The NGC 2024 proplyds (known so far) are all within 
0.2pc of the exciting UV sources. The mass loss rates of 
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Figure 6. Additional candidate proplyds in NGC 2024, which are fainter or less obvious than the main proplyds introduced here. Candidates 
5-7, detected in Paschen o are all close projected proximity to one another and pointing towards IRS 2b. Candidate 8 is detected in 
Ha and pointing towards IRS 1. These candidate proplyds are typically at least a factor 2 more distant from the UV sources than the 


unsubtle proplyds in the region. Note that the epoch is J2000. 


proplyds near IRS 2b are generally higher than those near 
IRS 1. This is at least in part due to the stronger UV flux 
from the higher mass IRS 2b. However the discs near IRS 
2b are also younger and therefore likely to be more massive 
and possibly also more extended which would also drive up 
the mass loss rate (e.g. Haworth et al. 2018a,b; Winter et al. 
2020b). 

'The mass loss rates are generally higher than those in- 
ferred for the NGC 1977 proplyds (~ 1078 Mo yr!) by Kim 
et al. (2016), with the exception of candidate proplyd 8. 
'This is because the NGC 2024 proplyds are closer to the UV 
sources (meaning stronger UV fluxes) as well as younger and 
hence possibly more massive/extended. 

'The ratio of a disc mass estimate to the mass loss rate 
gives a zeroth order estimate of the disc lifetime, however 
this comes with a number of caveats. The mass loss rate 
is typically a strong function of disc radius (e.g. Haworth 
et al. 2018b) and an evaporating disc gets truncated over 
time. High mass loss rates hence quickly result in a smaller 
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disc with a much lower mass loss rate that survives much 
longer than one would anticipate based on the earlier mass 
loss rate. In addition, the UV flux that any given disc is 
exposed to at < 0.1pc distances from the UV source for 
^kms ! velocity dispersions can change by up to two or- 
ders of magnitude within 1 Myr (Stórzer & Hollenbach 1999; 
Winter et al. 2019). We also estimated the mass loss rate us- 
ing the projected separation, when the true separation could 
be somewhat larger (and hence mass loss rate lower). The 
other key caveat here is that the van Terwisga et al. (2020) 
masses are in the mm dust only and it is the gas that is 
predominantly evaporated (Facchini et al. 2016; Sellek et al. 
2020). These caveats aside, assuming a dust-to-gas mass ra- 
tio of 1/100 the proplyd M/ M values are typically only a few 
tens of kyr, with the largest value being ~ 120 kyr for Pro- 
plyd 1. For the younger proplyds near IRS 2b the argument 
could be made that they have only recently begun being 
subject to photoevaporation (e.g. as IRS 2b ceased being so 
embedded along certain lines of sight). However proplyd 1, 
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Figure 7. The predicted ionisation front radius for a proplyd as 
a function of the separation fron IRS 1 (upper panel) and IRS 
2b (lower panel). The lines are our nominal values based on the 
ionising fluxes of Sternberg et al. (2003) and the shaded area 
goes out to ionising fluxes a factor 5 lower. In the upper panel 
the points are the 3 new proplyds in the vicinity of IRS 1 and 
the squares in the lower panel are new proplyds discovered in the 
vicinity of IRS 2b. 


near IRS 1, is in the ~ 1 Myr old region, raising the proplyd 
lifetime problem known from the ONC (e.g. O'dell & Wen 
1994). The solution to this seems to require either a lower 
dust-to-gas mass ratio (and hence higher gas mass than an- 
ticipated from continuum measurements, Sellek et al. 2020), 
ongoing star formation throughout NGC 2024 (even in the 
depleted western region) or that the remaining proplyds in 
the west recently migrated in towards IRS 1 (Winter et al. 
2019). 'The former solution of higher-than-canonical mass 
discs was also required by Haworth et al. (20182) to explain 
the formation of Trappist-1, even in weak UV environments 
and high mass discs have since been shown to be sustainable 
in a gravitationally stable manner around low mass stars 
(Haworth et al. 2020; Cadman et al. 2020). 


It is also important to note that proplyds have now 
been found in regions with ages ranging from « 0.5 Myr 
(this paper) to around 1.5 Myr (Kim et al. 2016; Kounkel 
et al. 2018), so higher gas masses, ongoing star formation, or 
continuous migration of discs into the high UV environment 
would have to be taking place to observe proplyds over such 
a length of time. 


6.2 Indirect evidence for external photoevaporation? 


Much recent evidence for disc evaporation has come indi- 
rectly from the distribution of disc properties in a cluster, 
for example searching for correlations between the disc mass, 
radius or disc fraction as a function of projected separation 
from a UV source (e.g. Mann et al. 2014; Guarcello et al. 
2016; Ansdell et al. 2017; Eisner et al. 2018). Haisch et al. 
(2000) found no obvious trend in the infrared disc fraction as 
a function of projected radius in NGC 2024, but we looked 
for any imprint of radiation environment on the van Ter- 
wisga et al. (2020) disc masses. 

NGC 2024 is a challenging region to look at in a statis- 
tical sense as proplyds are associated with both IRS 1 and 
IRS 2b. There is also a lot of gas and dust in the eastern 
region. The UV structure of the region is hence likely set 
by a complicated spatially varying combination of the two 
stars. 

Figure 8 shows the van Terwisga et al. (2020) disc dust 
mass estimates as a function of projected separation from 
IRS 1 and IRS 2b. The blue points represent the new pro- 
plyds introduced in this paper (where mass estimates are 
available) and the red arrows are the discs for which van 
Terwisga et al. (2020) could only place upper limits on the 
dust mass. We evaluate the Pearson's measure of linear cor- 
relation and Spearman's measure of higher order correla- 
tion. In each case we also evaluate the p-value by randomly 
shuffling the arrays of projected separation/disc mass ten 
thousand times, re-computing the correlation metrics and 
keeping track of the fraction of distributions with weaker 
metrics than the true one. 

In Figure 8 we include all discs with van Terwisga et al. 
(2020) mass constraints in each panel, regardless of whether 
or not they are closer in projected separation to IRS 1 or IRS 
2b. Disc masses are found to weakly increase with increasing 
distance from IRS 1 (Pearson's 0.2, p-value 0.12) (Spear- 
man's 0.27, p-value 0.05) though this will be influenced by 
the fact that there is a distinct population of younger more 
massive discs in the vicinity of IRS 2b. The disc masses are 
also found to weakly decrease with projected separation from 
IRS 2b (Pearson's -0.27, p-value 0.05), (Spearman's -0.2, p- 
value 0.13). The proplyds do not obviously stand out from 
the rest of the population, other than that they are some 
of the more massive members of the group (which gives a 
higher mass loss rate). 

In Figure 9 we again plot the disc mass as a function of 
projected separation, but filtering the data such that only 
the discs nearest to the UV source in projected separation 
are included. Correlation factors for IRS 1 are now (Pear- 
son's 0.028, p-value 0.81) (Spearman's 0.39, p-value 2x 1074) 
and for IRS 2b (Pearson's -0.23, p-value 0.013) (Spearman's 
-0.144, p-value 0.16). Most of the upper limits correspond 
to the same disc mass, but re-computing after discarding 
all upper limits less than 2.5 Mẹ only makes the Pearson’s 
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Figure 8. Disc mass as a function of projected separation from 
IRS 1 (upper panel) and IRS 2b (lower panel) for the discs with 
mass estimates from van Terwisga et al. (2020). The blue points 
are the proplyds introduced in this paper for which we could find 
van Terwisga et al. (2020) counterparts. The red arrows are the 
discs for which van Terwisga et al. (2020) have upper limits. 


and Spearman's correlations even weaker. So there is still 
no strong correlation with disc properties as a function of 
projected distance from either source. 

A key point to explain the lack of spatially varying en- 
vironmental imprint on the disc statistics is that the field 
considered here is rather small, with all discs being within 
~ 0.3pc of the UV sources. The correlation between disc 
mass and projected separation from ø Ori discovered by 
Ansdell et al. (2017) only kicks in beyond around 0.5 pc, 
with a flat distribution of masses interior to that. Over such 
a small distance scale the crossing timescale for a ~ 1 kms"! 
velocity dispersion is ^ 0.1 Myr. Since both the East and 
Western zones are many crossing timescales old the disc 
masses would be randomised with projected separation. 

In the older western region a larger fraction of discs only 
have upper limits on the mass, consistent with the picture 
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Figure 9. As in Figure 8, these plots show the disc dust mass as 
a function of projected separation from IRS 1 or IRS 2b. The 
difference here is that we only include the discs nearest to IRS 1 
(upper panel) or IRS 2b (lower panel) iin projected separation, 
whereas in Figure 8 all discs are included in both panels. The 
numbered blue points are the proplyds introduced in this paper 
for which we could find van Terwisga et al. (2020) counterparts. 
The red arrows are the discs for which van Terwisga et al. (2020) 
have upper limits. 


that this is older and more of the discs have been heavily 
depleted by photoevaporation. In the younger eastern region 
around IRS 2b, many discs are younger and are probably still 
well shielded by the immediate interstellar medium. It is also 
worth noting that all of the proplyds associated with IRS 2b 
are at the edge of the dense ridge, which would be consistent 
with the young stars/discs being progressively exposed to 
the UV as the gas that they are embedded in is dispersed. So 
even if the region were not old enough for the crossing time 
to wash out any spatially varying imprint of evaporation, 
there may be no signature of evaporation as a function of 
projected separation simply because the discs are young and 
embedded. 
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6.3 External photoevaporation in NGC 2024: reflecting 
upon the van Terwisga et al. (2020) interpretation 


van Terwisga et al. (2020) presented ALMA 1.3mm contin- 
uum observations of the two distinct populations in NGC 
2024. They suggested that a key distinction was that the 
mm-depleted western population is much less embedded and 
possibly subject to evaporation by IRS 1 and IRS 2b. They 
also suggested that the discs in the western population may 
have just formed smaller and hence had shorter viscous life- 
times. It is important to understand if it is possible for there 
to be systematic regional variations in disc properties in- 
dependent of environmental influence like photoevaporation 
(e.g. set by the initial molecular cloud properties there). One 
possible mechanism to enable this would be variations in 
the cosmic ray ionisation rate, which Kuffmeier et al. (2020) 
showed can affect the magnetorotational braking during col- 
lapse and lead to smaller disc sizes. 

Our discovery of proplyds in NGC 2024 has some im- 
portant implications for the discussion above. First of all, 
we have indeed confirmed that external photoevaporation is 
at work in the region. In the western region our proplyds 
point towards IRS 1, implying that the IRS 1 dominates the 
UV radiation field there and could have depleted the disc 
fraction in the region even without the help of IRS 2b. 

What was not anticipated from van Terwisga et al. 
(2020) is ongoing photoevaporation in the eastern popula- 
tion, but with 5 of our proplyds in this zone we have shown 
that this is actually the case. Despite the eastern region be- 
ing younger and having a higher disc fraction in terms of 
ALMA detections (45 +7 per cent) compared to the western 
region, the disc fraction is still low compared mm contin- 
uum detection rate in other star forming regions, such as 
Lupus (70per cent mm continuum detection fraction in a 
1-2 Myr region Ansdell et al. 2016) and Sigma Ori (40 per 
cent mm continuum detection fraction in a 3-5 Myr region 
Ansdell et al. 2017). Though future comparisons ensuring 
equivalent detection thresholds will be necessary, this is ev- 
idence that many discs in the eastern region of NGC 2024 
that have not been shielded by the ISM have been subject 
to photoevaporation by IRS 2b. 


6.4 Implications for planet formation 


Our understanding of planet formation in the idealised sce- 
nario of discs as isolated systems is still not yet fully under- 
stood. By extension, the added complication of the impact 
of external photoevaporation on planet formation is not well 
understood. What is known is that external photoevapora- 
tion reduces the disc mass (e.g. Ansdell et al. 2017; Winter 
et al. 2018a; Sellek et al. 2020), radius (e.g Haworth et al. 
2017; Eisner et al. 2018; Sellek et al. 2020) and lifetime (e.g. 
Guarcello et al. 2016; Concha-Ramírez et al. 2019b; Winter 
et al. 2020a) which could all of course indirectly affect the 
formation of planets. 

A key element for planet formation, even in the inner 
disc, is the dust reservoir of the entire disc. In the outer 
disc grains grow and drift inwards to smaller radii and can 
contribute to planet formation there, which is expected to 
be necessary to form the large mass of close-in planets (e.g. 
Birnstiel et al. 2012; Marcy et al. 2014; Ormel et al. 2017). 

In the outer disc, small dust grains get entrained in the 


photoevaporative wind (Facchini et al. 2016), whereas larger 
grains are pushed inwards with the retreating disc outer edge 
and could still ultimately contribute to planet formation 
(Sellek et al. 2020). The total dust mass lost in a disc is set by 
the competition between grain growth in the disc (which is 
fastest at small orbital radii) and the rate of external photo- 
evaporation (e.g. Haworth et al. 2018a; Sellek et al. 2020). If 
grain growth proceeds through the disc quickly enough, then 
only a small amount of dust would end up entrained in the 
wind. Conversely very strong, early photoevaporation could 
deplete significant amounts of the potential dust reservoir 
that could otherwise have migrated in to contribute to inner 
planet formation. A key question is therefore how long discs 
remain embedded, permitting dust to grow to a size where 
it is immune to photoevaporation once the source becomes 
exposed to the ultraviolet source. 

In light of the above, our finding evidence for very early 
photoevaporation throughout NGC 2024 has added impor- 
tance because it may also be significantly reducing the dust 
mass reservoir (in addition to the usual expected impact on 
the gas disc mass/radius/lifetime). This would also explain 
why the dust mass estimates of discs in the 0.2 — 0.5 Myr 
eastern population of NGC 2024 are more comparable to 
those in 1 — 3 Myr old regions (van Terwisga et al. 2020). 

There is growing evidence for planet formation on < 
1 Myr timescales (e.g. Segura-Cox et al. 2020), but here we 
show evidence for very early external photoevaporation (0.2- 
0.5 Myr). 


6.5 Potential for additional evaporating discs 


Our search for externally photoevaporating discs in NGC 
2024 using archival data is certainly not complete. The 
F187/F190 fields using NICMOS have a small field of view 
and the archival data does not cover the whole region. Future 
HST WFC3 Ha and Paschena surveys may reveal more 
proplyds in the central region. Additional searches at larger 
distances from the less embedded western region may also 
reveal photoevaporative winds from discs where the UV field 
is weaker, e.g. in CI (Haworth & Owen 2020). Surveys out 
to larger distances from IRS 1 and IRS 2b may also re- 
veal any statistical imprint of photoevaporation on the disc 
properties. The above combined with the need to confirm 
our candidate proplyds means that NGC 2024 is a prime 
region to search for further evidence for external disc pho- 
toevaporation. 


7 SUMMARY AND CONCLUSIONS 


NGC 2024 is a young star forming region in Orion with 
a radial age gradient (Getman et al. 2014). van Terwisga 
et al. (2020) also recently found two slightly differently 
aged sub-populations with very different ALMA 1.3mm 
continuum disc fractions in the inner (« 0.5pc) part of 
the region. To gauge the role of external photoevaporation 
of discs in sculpting these populations we searched for 
proplyds in archival HST observations of NGC 2024. We 
draw the following main conclusions in this work. 


1) We discovered 4 firm proplyds and 4 candidate proplyds. 
Three of these are in the vicinity of the B0.5V star IRS 
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1 and the other 5 associated with the (younger and more 
embedded) O8V star IRS 2b. External photoevaporation of 
discs is definitely occurring in NGC 2024 and within both the 
statistically older and younger populations. It is also worthy 
of note that together with NGC 1977 (Kim et al. 2016) 
early type B stars have now clearly been demonstrated as 
capable of evaporating the discs of nearby stars. 


2) The geometry of the cometary proplyds introduced here 
can be explained using only IRS1 and IRS2b. No further 
embedded UV source is required. Proplyds 1 and 2 in the 
western region show evidence for a jet because of elongated 
emission through the cusp. 


3) Proplyd 4 in this paper is possibly a binary proplyd with 
a separation of at least 400 au, though we cannot discern an 
interproplyd shell with the archival Paschen a observations. 
Follow up observations of this system will be particularly 
interesting. 


4) With a maximum proplyd M/M of ~ 120kyr, even 
with the caveats of such an estimate there appears to be a 
proplyd lifetime problem in NGC 2024. Furthermore, the 
proplyds are actually some of the more massive discs from 
van Terwisga et al. (2020). This suggests that either discs 
are more massive in gas than expected from continuum 
measurements (ie. that the dust-to-gas mass ratio is 
lower than 1/100 Sellek et al. 2020), that there is ongoing 
migration of relatively massive discs into the high UV 
inner part of the cluster, or ongoing star formation, in both 
regions (Winter et al. 2019). 


5) There is no clear statistical imprint of evaporation on 
the disc masses in the region as a function of projected sep- 
aration from the main UV sources, but this is unsurprising 
given the young age of the region and that the field being 
considered is very small so discs that are not embedded 
are evaporated to the same degree (i.e. rapidly). Further- 
more, the dynamical time in this compact central region 
is short (~ 0.1 Myr) so any imprint of photoevaporation 
with projected separation would also quickly be washed out. 


6) Our results slightly change the van Terwisga et al. 
(2020) interpretation of the region in a few ways. Firstly, 
evaporation of the depleted western population appears to 
be dominated by IRS 1, since the proplyds in the region 
are pointing towards it, and so IRS 2b is not required to 
play a role there. Second, the eastern population is also 
being actively evaporated by IRS 2b where discs are not 
embedded in the ISM, which may also explain why the 
ALMA mm continuum disc fraction is more like that of an 
older region (van Terwisga et al. 2020). 


7) Our evidence for external photoevaporation in the 
0.2-0.5 Myr eastern region shows that this environmental 
influence can be in competition even with very early planet 
formation (e.g. Segura-Cox et al. 2020). 
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